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As part of our ongoing Wolf-Rayet (WR) Magellanic Cloud survey, we have discovered 13 new 
WRs. However, the most exciting outcome of our survey is not the number of new WRs, but 
their unique characteristics. Eight of our discoveries appear to belong to an entirely new class 
of WRs. While one might naively classify these stars as WN3+03V binaries, such a pairing 
is unlikely. Preliminary CMFGEN modeling suggests physical parameters similar to early-type 
WNs in the Large Magellanic Cloud except with mass-loss rates three to five times lower and 
slightly higher temperatures. The evolution status of these stars remains an open question. 


1 Discovery 

We recently began a multi-year project surveying 
the Magellanic Clouds for Wolf-Rayet (WR) stars. 
Using the Swope 1-m at Las Campanas, we took im¬ 
ages using on- and off-band interference filters and 
then used a image subtraction program to identify 
WR candidates. The candidates were then spectro¬ 
scopically confirmed using MagE on Magellan. So 
far we have covered ~60% of the Magellanic Clouds 
and discovered 13 new WRs (Massey et al. 2014, 
2015). But, it isn’t the number of newly discovered 
WRs that makes our survey so successful, but the 
types of these stars. Eight of them appear to belong 
to an entirely new class of WRs. 

Based on the spectrum shown in Figure 1, one 
might naively classify these stars as WN3+03V bi¬ 
naries. The WN3 classification comes from the star’s 
Nv emission (AA4603,19 and A4945), but lack of 
N IV. The 03V classification comes from the strong 
Hell absorption lines but lack of Hel. However 
such a pairing of an WN3+03V is highly unlikely. 
First, 03Vs are quite rare since they are the hottest 
and most luminous of the dwarfs (Ekstrom et al. 
2012). In the Large Magellanic Cloud (LMC), only 
around a dozen are known outside of the 30 Dor 
region (Skiff 2014). Second, assuming single star 
evolution, the lifetimes of these stars preclude such 
a pairing: WN3s take ^3-5 million years to form 
whereas a massive star evolves out of the 03V stage 
after only a million years. Third, we have yet to see 
any radial velocity variations as you would expect 
in a binary system, though we are still gathering 
data for confirmation. Fourth, our UV data (dis¬ 
cussed later) shows no C IV A1550, the most promi¬ 
nent UV line in O-stars. However, the most convinc¬ 
ing evidence deals with the stars’ absolute magni¬ 
tudes. The WN3/03Vs we have observed are faint 
with My ~ —2.5. 03Vs, on the other hand, have 
much brighter absolute magnitudes with My ^ — 5.5 
(Conti 1988). Thus, we conclude these stars are not 
WN3+03V binaries and adopt the naming conven¬ 
tion of WN3/03V. 
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Fig. 1: Spectrum of LMC170-2, one of our newly dis¬ 
covered WN3/03V stars (Massey et al. 2015). 


In addition to the MagE spectra (resolving power, 
R , of 4100) obtained for all eight newly discovered 
stars, we observed three stars in the UV using COS 
on HST. We also observed one star using both FIRE 
and MIKE on Magellan, giving us the NIR and an 
optical spectrum with R = 11000 — 14000. Next we 
attempted to model our spectra using a single set of 
physical parameters. 


2 Physical Parameterization 

To model the WN3/03s, we turned to CMFGEN 
(Hillier & Miller 1998), a stellar atmosphere code 
designed for hot stars with stellar winds where the 
usual assumptions of plane-parallel geometry and 
LTE no longer hold. Using only the MagE optical 
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data, we produced a good fit to both the emission 
and absorption lines in the spectra, as shown in Fig¬ 
ure 2. Table 1 lists the physical parameters used. 



Fig. 2: Best CMFGEN fit to the optical spectrum of 
LMC170-2 from Massey et al. (2014). 


Tab. 1: CMFGEN model parameters for LMC170-2 


Parameter Best Fit 


T eff (K) 

L/L e 

M* (M 0 /yr) 
He/H (by #) 
N 

C, O 


100,000 
4 x 10 5 
1.2 x KT 6 
1.0 

10.0 x solar 
0.05 x solar 


Note: Assumes a clumping filling factor of 0 . 1 , Uoo =2400 km s 1 , 
/ 3 = 0 . 8 , and u s i n i = 150 km s -1 . 

The majority of these parameter values are com¬ 
parable to those found for early-type LMC WNs 
(Hainich et al. 2014). Despite the stars’ faint vi¬ 
sual magnitudes, their bolometric luminosities are 
normal. These stars appear to be evolved, with sig¬ 
nificantly enriched N and He. While their effective 
temperatures are a bit on the high side, they are still 
within the range found for other WN3s. The most 
unusual value is their mass-loss rate. As Figure 3 
shows, our WN3/03 star LMC170-2 has a mass- 
loss rate that is three to five times lower than those 
found for other early-type LMC WNs as analyzed 
by Hainich et al. (2014). As discussed later, this 
low mass-loss rate appears to hold true for all eight 
WN3/03s currently known. 
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Fig. 3: Luminosity vs. mass-loss rate for LMC170-2 
(large red star) compared to those of other early-type 
LMC WNs analyzed by Hainich et al. (2014). Adapted 
from Fig. 6 of Hainich et al. (2014) and used with per¬ 
mission. 


Since the physical parameters for this model were 
determined using only the optical data, we were cu¬ 
rious to see how well the model would match our 
newly obtained UV and NIR spectra. As is shown in 
Figure 4, the agreement between the model and our 
UV and NIR spectra is practically perfect. Due to 
space constraints we have only shown the model fit 
for He IIA1640 in the UV, but it should be noted that 
all of our diagnostic lines such as the OIVA1038 reso¬ 
nance doublet, N VA1240, and the lack of C IVA1550 
are similarly well fit. As an added bonus, the spec¬ 
tral energy distribution of the model almost per¬ 
fectly matches our UV, optical and NIR spectra. 
This confirms that the reddening is well determined. 
However, it should be noted that for stars with 
T e ff > 30,000 K, the spectral energy distributions 
will all look similar longwards of 1000A as this is 
the Rayleigh-Jeans tail of the black body distribu¬ 
tion. 
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Fig. 4: Best CMFGEN fit to the UV and NIR spectrum 
of LMC170-2. 
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With our exemplary model in hand, we began a 
sensitivity analysis of the model’s input parameters. 
We first attempted to slightly raise and lower the 
value of each parameter (e.g., increase by 10%, then 
decrease by 10%). As co-author John Hillier wrote, 
we soon discovered that, “it’s an interesting regime 
- M, H/He, T e ff, and log g all seem to interact in 
complex ways.” For example, decreasing T e ff from 
100,000 K to 90,000 K yielded a clearly inferior fit 
to the He and H lines. However, increasing T e ff to 
110,000 K brought the model too close to the Ed¬ 
dington limit. This could be avoided by changing 
the surface gravity in parallel with T e ^. So, doing a 
clear-cut sensitivity analysis proved to be more dif¬ 
ficult than originally expected. Still, we were able 
to get reasonable uncertainty estimates for the ma¬ 
jority of the parameters. However, we hope to refine 
these once we have separately modeled all eight of 
the spectra. 

Figure 5 shows the striking similarities between 
the spectra of our WN3/03s, suggesting that a good 
model fit to one star will also fit the remaining seven 
stars quite nicely. And indeed, this is what we found. 
While there are slight differences, we expect that 
the physical parameters will only change by small 
amounts for each star. Certainly none of the spectra 
suggest a drastically different temperature or mass- 
loss rate (for example). 



Fig. 5: Spectra of five of our newly discovered WN3/03s 
from Massey et al. (2014). 


3 Next Up 

We still have a lot to learn about these objects and 
certainly there are many questions left to be an¬ 


swered. We are currently working on constraining 
the model parameters and determining whether the 
“best fit” model presented above is unique or if a 
similarly good model can be obtained using a dif¬ 
ferent set of parameters. Once this has been deter¬ 
mined, we will continue modeling all eight stars to 
get a more global view of the range of expected pa¬ 
rameter values within this WR subtype. At that 
point we can begin working with our stellar evo¬ 
lutionary colleagues to determine whether current 
models predict the existence of these stars. 

We additionally hope to get a better handle on 
where these stars come from. If they were spa¬ 
tially separated, or more isolated, from the “nor¬ 
mal” LMC WRs, one might infer that the WN3/03s 
were formed from different (lower metallicity?) pro¬ 
genitors. However their spatial distribution suggests 
that this is not the case. We also haven’t ruled out 
an explanation dealing with binary evolution. While 
they don’t appear to be binaries now, it is possible 
that they were in the past. As we continue on with 
our WR survey in the Magellanic Clouds, we hope to 
find even more of these exciting new stars. Within 
the next year we should have a much better grasp 
of their numbers, physical parameters and possibly 
even their origins. 
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